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ABSTRACT 

Over a broad range of initial inclinations and eccentricities an appreciable fraction of hierarchical triple 
star systems with similar masses are essentially unaffected by the Kozai-Lidov mechanism (KM) until the 
primary in the central binary evolves into a compact object. Once it does, it may be much less massive than 
the other components in the ternary, enabling the "eccentric Kozai mechanism (EKM):" the mutual inclination 
between the inner and outer binary can flip signs driving the inner binary to very high eccentricity, leading to 
a close binary or collision. We demonstrate this "Mass-loss Induced Eccentric Kozai" (MIEK) mechanism by 
considering an example system and defining an ad-hoc minimal separation between the inner two members at 
which tidal affects become important. For fixed initial masses and semi-major axes, but uniform distributions of 
eccentricity and cosine of the mutual inclination, ~ 10% of systems interact tidally or collide while the primary 
is on the MS due to the KM or EKM. Those affected by the EKM are not captured by earlier quadrupole-order 
secular calculations. We show that fully ~ 30% of systems interact tidally or collide for the first time as the 
primary swells to AU scales, mostly as a result of the KM. Finally, ~ 2% of systems interact tidally or collide 
for the first time after the primary sheds most of its mass and becomes a WD, mostly as a result of the MIEK 
mechanism. These findings motivate a more detailed study of mass-loss in triple systems and the formation of 
close NS/WD-MS and NS/WD-NS/WD binaries without an initial common envelope phase. 
Subject headings: stars: binaries: close, — celestial mechanics, stellar dynamics — stars:white dwarfs — 
supernovae: general 
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1. INTRODUCTION 

iKozail d 19621) and lLidovl d 19621) showed that a hierarchical 
triple system with an inner binary of masses mo and mi, and 
a tertiary of mass ni2, can exchange angular momentum be- 
tween the inner and outer orbits periodically. These oscilla- 
tions drive the eccentricity of the inner binary (ei) to high 
values if the tertiary is inclined with 39.2° < i < 141.8°. 
For nearly circular orbits, in the test-particle approxima- 
tion (mi <^ mo,m2), the maximum eccentricity is given by 
(llnnanenet al.ll 19971) 



in. max — I 1 ^ COS i 



1/2 



(1) 



when the three-body Hamiltonian is expanded to quadmpole 
order in the ratio of the inner to outer semi-major a xes (ai/a^X 
and neglecting tidal forces and general relativity (|Blaes et aTl 
l2002tlMiller & Hamiltonll2002tlFabrvckv & Tremaindl2007l) . 

Tidal friction in Kozai-affected triples can cause a\ to 
decrease. As the inner binary is driven to high ec- 
centricity, the periastron of the secondary approaches the 
primary, and tidal frict ion tends to circularize the orbit 
dMazeh & Shaharr] fl979). This mechanism has been pro- 
posed as an explanation for the prevalence of "hot" Jupiters 
with few-day orbits around their host stars (|Wu & Murrav 



2003HFabrvckv & Tremainel2007tlWu et al.l2007l:lNaoz et alj 
2011 alibi) , the very high triple fraction of close solar-type bi- 
naries with periods less than ~ 5 days dTokovinin et al.ll2006t 
Fabrv ckv & Tre maine 20071 ), and the formation of blu e strag- 
glers in globular clusters dPerets & Fabrycky 2009). The 
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combination of Kozai cycles and tidal friction has also been 
explor ed ge nerally for triple star systems i n iKiseleva et al.1 
(1998) and lEggleton & Kiseleya-Eggletonl d2001l) . As dis- 
cussed by IWu & Murrayl (120031) . for the case of Jupiter-mass 
planets around Sun-like stars, by equating the tidal and Kozai 
forcing terms in the equation for the secular evolution of e\, 
one estimates the circularization semi-major axis of the inner 
bi nary to be ~ 3R(ry 

ILithwick & Naozl (1201 ll) . iKatz et all d2011l) . and 
iNaoz et all d201 lalfa) have recently emphasized that the 
commonly used quadrupole-order expansion of the three- 
body Hamiltonian is insufficient to capture the secular 
dynamics of triple systems in the test particle approximation 
when the outer eccentricity (e 2 ) is non-zero. The parameter 



/ mo-mA ( aA e 2 
\ mo + mi ) \a 2 ) 1-e 



(2) 



measures the importance of the octupole-order terms in 
the doubly-averaged thr ee-body Hamiltoni an relative to the 
quadrupole-order terms. INaoz et ail (1201 lab first showed that 
when (mi <C mo,m 2 ) and e 2 ^ 0, it is possible for the triple 
to "flip": the system exhibits quasi-periodic cycles in cos i 
through 0, and the tertiary passes from prograde to retrograde 
and vice versa. These flips occur even for e oct as small as 
10~ 3 in the test particle approximation, provided that the sys- 
tem is sufficiently inclined initially and the arguments of pe- 
riastron of the inner orbit (g\) and the o uter orbit (#2) are 
chose n ju diciously (see F igures 7 and 3 of ILithwick & Naozl 
1201 ll and lKatz et aill201 ll respectively). During these flips, 
the inner binary is driven to extremely high eccentricities 
(\-e\ ~ 10~ 5 ). This behavior occurs for a broad range of 
ohlo parameiers in octupole-order calculations, but it is not present 
in the quadrupole-order calculations, and has been referred to 
as the "eccentric Kozai mechanism." 
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An important piece of physics missing from the test particle 
approximation of t he eccentric Kozai mechanism i s the (mp — 
mi) t e rm in e oct ( Krvmolowski & Mazehl fl 999: Fo rd et all 
I2000t iBlaes et alj 120021; iNaoz et all 1201 lbl) . Because the 
three-body Hamiltonian in Jacobi coordinates, expanded to 

all orders in 02/01 has a sum over (m J - (-mi) 7 " 1 ) for j = 
2,3, • • • , when mo = mi all j = odd terms are zero (Harrington 
1968). Importantly, the flip phenomenon present in the test- 
particle calculations at modest ei and i is strongly suppressed 
as mo — > ni\ because the octupole-order terms become neg- 
ligible 3 . Our own experim ents with the octupole -order code 
based on lBlaes et all J2002I) and lThompsonl ( 1201 ll) shows that 
in order for flips to occur over a broad range in g\ and g2, mi 
must be less than mo by a factor of ~ 2 for a given 01/02 and 

The dependence of the flip phenomenon on (mo -mi) im- 
plies that many systems will be affected by mass-loss as the 
primary becomes a WD or a NS. Indeed, many MS triple 
star systems will have mo ~ m\, because the inner bina- 
ries of solar-like triples typically have a flat mass distribu- 
tion dMazeh et alJl!992l : lRaghayan e t al. 20 1 0]) and may even 
have a preference for "twins" (Pinsonneault & Staneki l2006l 
Ragha van et al.|[20ToT) . The similar masses of the inner binary 
suppresses the eccentric Kozai mechanism since e oct is small. 

Furthermore, a significant fraction of triple systems (those 
with modest mutual inclination) will have inner binaries that 
are wide enough to avoid tidal contact and collisions on the 
MS, even if the system undergoes normal Kozai-Lidov oscil- 
lations with <? max given approximately by Equation £[). How- 
ever, once the primary of the inner binary evolves off the MS 
and loses the majority of its mass, the system effectively en- 
ters the test-particle approximation in many cases, particularly 
for intermediate- and high-mass primaries where the change 
in mass as the MS star becomes a compact object is large. 
After mass loss, e oct can increase dramatically, allowing the 
systems to exhibit the eccentric Kozai mechanism, and driv- 
ing the inner binary to extremely high eccentricities and tidal 
contact or collision. The likely result of the MIEK mecha- 
nism is a close compact-object(WD/NS)-MS binary with sep- 
aration ~ R . Such a configuration would produce a variety 
interesting astrophysical systems. 

In this paper we demonstrate this "Mass-Loss Induced Ec- 
centric Kozai" (MIEK) mechanism and provide a preliminary 
exploration of its dependence on the initial eccentricities and 
mutual inclination of the triple system, saving a detailed treat- 
ment for a future paper. In Section |2] we describe our method 
for integrating the orbits of triple systems with mass-loss, and 
we give two example systems with different mass-loss time- 
scales to illustrate the MIEK mechanism. By directly inte- 
grating triple systems with an Af-body code we circumvent 
the need to pick a limiting order when expanding the three- 
body Hamiltonian and solving the secular dynamics. In Sec- 
tion|4] we explore parameter space by investigating the effects 
of eccentricity and inclination on the MIEK mechanism for 
an example ternary. We find that a fraction of widely sepa- 
rated triple systems are brought to tidal contact for the first 
time only after the primary evolves off the MS. In Section 
[5] we discuss the possible outcomes and the implications of 
systems whose inner binary comes to tidal contact or collides 

3 Intuitively, any odd-order terms in the multipole expansion of the three- 
body Hamiltonian are odd functions of roi, the separation vector between 
members of the inner binary, and thus averages to zero for symmetric inner 
masses. 




FIG. 1 . — Comparison between our modified FEWBODY code (black), and 
an octupole-order secular code (red, Blaes et al. 2002). The triple system has 
m = 7.OM ,m, =6.5 M , m 2 = 6 M Q , a x = 10 AU, a 2 = 250 AU, e x =0.3, 
«2 - 0.1, gi = 0°, g2 = 0°, and cos i = 0.4. There is a good agreement between 
the calculations but there is a slow drift in gi —gj. Figure discussed in i|2] 

at each evolutionary phase. In Section [6] we provide a brief 
conclusion. 

2. SIMULATING MASS-LOSS IN TRIPLES 

To simulate the dynamical effects of mass-loss on hierarchi- 
cal triple systems, w e modified the A^-body code FEWBODY 4 
(Frege auet al.l20 04) to create triple systems and evolve them. 
We create a triples with specified inner and outer semi-major 
axess (ojand 02), eccentricities (<?iand <?2), arguments of pe- 
riastron (giand gi), primary mass in the central binary (mo), 
secondary mass in the central binary (mi), tertiary mass (m2), 
random phase angles, and mutual inclination between the two 
orbits (cos 2). 

During integration of the triple, we check that the triple 
remains bound and stable. If the triple becomes unbound 
we continue integrating the system until a conservative min- 
imum tidal perturba tion on the binary is met (S < 10~ 7 ; see 
Frege auet al.l l2~004). and we thus include any possible reso- 
nant interactions which could result when the triple system 
becomes unstable 5 . Integration is also stopped if the time 
steps become exceedingly small. This occurs for a handful 
of triples where the separation between the inner binary ap- 
proaches because the stars are on nearly radial orbits. How- 
ever, in the current study these cases only occur when the in- 
ner two stars are separated by -C Rq . 

To check our modified FEWBODY integrator we tested a 
sample of triple systems against the N -body integrator MER- 
CURY ( Chambers & Migliori nil 19971) and found good agree- 
m ent. We then tested our code against the octupole-order code 
of IBlaes eT al. (2002) with the general relativistic terms re- 
moved. An example triple system is shown in Figure Q] As 
can be seen in the figure there is reasonable agreement with 
the secular octupole-order calculation. There is, however, a 
slow drift in g\ —gi between the N-body codes and the secular 
code. 



4 FEWBODY is now available at http://fewbody.sourceforge.net/ 

5 This allows our calculations to include systems affected by the "triple 
evolution dynamical instability" (TEDI.Perets & Krattei 20 3). 
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2.1. Evolutionary Phases 

When we evolve a triple system, we break its evolutionary 
phases into four distinct stages: (1) primary on the MS, (2) 
primary evolved into a giant but prior to significant mass loss, 
(3) during a period of significant mass loss, and (4) after the 
end of mass loss when the primary has become a compact ob- 
ject. However, because of the prohibitive computational time 
required to integrate close and highly eccentric triples for the 
full MS lifetime of the primary, we truncate o ur integrations 
after 10 3 Kozai-Lidov tim es (?k) defined by (llnnanen et alj 
[l997t lHolmaneFaII[T997h : 



4 / a\(mQ + m\) 



Gm\ 



;6.4xl0y(^fY^^Y7^) (-^ 
y UOAU/ V 13 - 5M o/ V m 2 J \25ai 



f MS = min(f ms , o, 10 3 f K ) x (0.95 + 0. 10, 



f RG = min(0.1f ms .o,10 t K ). 



m f = (0. 109 ± 0.007) m + (0.394 ± 0.025) M e 



?after = min(f„ 



l.lf/r 



-f m i, 10 ?k), 



(3) 



where b2 = ai(\—ei) x l' 1 - Thus, we implicitly assume that the 
triple systems have no interesting behavior after this time. 

We account for the MS lifetime of the primary by only in- 
tegrating the triple systems for a time (?ms) defined by: 



(4) 



where £ is a uniformly distributed random number on the 
interval [0,1), f mSj o is the MS lifetime of the primary. We 
determine the MS lifetimes by interpolating the logarith- 
mic lifetimes on a grid of stellar models run with the Yale 
Rotating Evolution Code with i nput physics described in 
Ivan Saders & Pinsonneaultl (1201 ll) . We multiply by the sec- 
ond factor in Equation (|4]i in order to randomize exactly when 
the primary evolves off the main sequence to avoid any cor- 
relations of the start of mass loss between different triple sys- 
tems in our sample. 

To simulate the giant phase of the primary, we continue to 
integrate the triple for a time (?rg) defined by: 



(5) 



Then, to simulate mass-loss, we linearly decrease mo over a 
mass-loss time scale (f m i) to the WD mas s as sp ecified by the 
initial-final mass relation of iKalirai et alj Q2008) given by: 



(6) 



where mo,f is final mass of the primary. To test our mass-loss 
prescription, we integrated binaries where one member under- 
went mass loss and tested that the change in the semi-major 
axis and eccentricity of the binary matched that for instanta- 
neous mass loss and adiab atic mass loss (i.e., Equations 4-6 in 
iBaribault & Pineault|[l994l) when f m i < P and when t m \ > P, 
respectively, where P is the orbital period of the orbit. We 
found good agreement. 

After mass-loss has finished we integrate the evolved triple 
for a time (f a fter) given by: 



(7) 



where f ms> i is the MS lifetime of the next most massive star 
in the triple system and t^ is redefined by Equation (O with 
the orbital parameters of the triple right after the conclusion 
of mass-loss. 

2.2. Tidal Contact and/or Collisions 



A detailed treatment of stellar tides and tidal dissipation is 
beyond the scope of this paper. However, to determine when 
tidal effects are likely to be important, we define an ad-hoc 
minimal separation between the inner two binary members, 
'"tide, below which strong tidal contact or a collision is assumed 
to have occurred. That is, if the separation between the pri- 
mary and secondary is r < r t a e , we assume either (1) that tidal 
dissipation will remove energy from the inner orbit decreasing 
a i or (2) if the separation between the inner binary members 
changes dramatically in a single inner orbital period (Pi), a 
physical collision may o ccur before tides bec ome important, 
as in the recent work of iKatz & Pond (120121) . In either case, 
our simulations are not valid for r < r t id e - 

When choosing a tidal criterion we considered two tidal 
effects which will alter the evolution of the system. The 
first is the nondissipative contribution that tidal bulges on the 
members of the inner binary have on apsidal motion (tSternel 
119391) . This addi tional apsidal motion may detune the Kozai 
mech anism (e.g. IW u & Murray 2003; Fabrycky & Tremainel 
120071) . We find that the timescale of this apsidal motion 
becomes comparable to t& when r per j ~ 4/? mSj o, assuming 
the typical tidal Love number valid for n = 3 polytropes, 
k = 0.028, (Eggleton & Kiseleva-Eggletori l2001l) and the sys- 
tem parameters for the fiducial system presented in §|3] The 
second effect is eccentric i ty damping due to dissipative tides 
dMazeh & Shahaml [19791 IWu et alj [20071) . For our fiducial 
system, we derive a critical inner seperation by equating the 
tidal and Kozai forcing terms in the equation for the secular 
evolution of e\. We find a critial inner separation of: 

/10 4 k \ 2/ ' 3 

rcnt ~ 3R ™<°{-Qom8) ' (8) 

where Q is the tidal dissipation factor. As can be seen in 
Equation ||8), r cr j t is very insensitive to the exact values cho- 
sen for k and Q. For simplicity, and to allow a broad, al- 
beit incomplete, survey of parameter space, we take r t id e = 
3 x max(R ms> o,Rms, l) for the bulk of this study, bu t we inves- 
tigate the effect of changing r t id e on our results in §14.31 

2.3. General Relativistic Precession 

General relativistic (GR) periastron precession can "de- 
tune" the normal Kozai mechanism, and truncat e the maxi- 
mum eccentricity attainable for a triple system (|Blaes et all 
12001 iMiller & Haifmtonl l200l iThompsonl 120111) . FEW- 
BODY does not contain post-Newtonian corrections, so 
we must check if it is permissible to ignore GR preces- 
sion. The GR precessio n timescale (e.g., Equation 23 of 
Fabrycky & Tremaine 2007) is given by: 



1 «i 

f GRp = — — 



a\c 



3 c \G(mo + m\) 



3/2 



:3.4xlOV 

\(m +mi) 



3/2 



fll \ 5 /2 

lOAu) 



(9) 

As discussed in Blae s et aT] d2002[) . the Kozai mechanism only 
operates if ta < ?GRp- We see that for the masses and semi- 
major axes of our fiducial case presented in ^3] ?k < ?GRp for 
any choice of <?2 until the inner eccentricity is increased such 
that r pel -i < R m ^ \ . Thus the normal Kozai mechanism is free 
to operate in our triple integrations. 

However, we also see in §f5] that the eccentric Kozai 
timescale (?ekX the timescale for "flips" and extreme eccen- 
tricity maxima, is much longer then ?k, ?ek ~ 10 6 years. 
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Furthermore, we see that ?ek > ?GRp for values of e\ where 
'"pen > Rq, which would seem to imply that GR precession 
might detune the eccentric Kozai mechanism before tidal ef- 
fects become important. This is not clear however, because 
the majority of the GR precession occurs when e\ is large, 
which is only a small fraction of the eccentric Kozai cycle. To 
test if GR precession will suppress the maximum e\ reached, 
we take the fiducial triple system after mass loss that exhibits 
the eccentric Kozai flip, and compare this to the equivalent 
octupole-or der calculat i on tha t includes GR prec ession using 
the code of lThompsonl d2011l) dBlaes et alj|2002h . We found 
good agreement between the two calculations, implying that 
GR precession does not shut off the eccentric Kozai mecha- 
nism, and is thus not important for the current study. How- 
ever, these effects should be more carefully considered in fu- 
ture studies where a\ is much smaller, and when the perias- 
tron distance can be significantly smaller than Rq if both the 
primary and secondary are compact objects. 

3. MIEK: 

THE MASS-LOSS INDUCED ECCENTRIC KOZAI MECHANISM 

Figure [2] presents two example systems with the left and 
right panels having mass-loss time scales of ? m i = 10 4 yr and 
t m \ = 10 6 yr, respectively. Both have mo = 7.0 M©, m.\ = 
6.5 M Q , m 2 = 6 M Q , a\ = 10 AU, a 2 = 250 AU, e x = 0.1, 
e 2 = 0.7, gi = 0°, g 2 = 180°, and cos i = 0.5 (z = 60°). The top 
panels present the evolution of a\, a 2 , e 2 , and mo. The mid- 
dle and bottom panels show the evolution of cos z and 1 — e\, 
respectively. In the bottom panel, the red line shows where e\ 
is high enough such that the radius of periastron of the inner 
binary (r per j) is equal to R©. The example systems go through 
normal Kozai-Lidov cycles before mass-loss begins, reaching 
a maximum e\ ~ 0.75 as predicted by Equation (fl}, but r per j 
is still more than two orders of magnitude too large for tidal 
interactions to be significant. 

After adiabatic mass-loss mo — > 1.15M Q , and e oct increases 
(Equation To calculate this increase analytically, we first 
note that the increase in the semi-major axis and change in 
eccentricity of a binary that underwent adiabatic mass-loss is 
given by: 

Of M Q 



— = — , and <?f : 

ao Mf 



' <?o, 



(10) 



where Mo and Mf are the initial total mass and final total mass 
of the binary, respectively. From Equation (TfOt we see that 
the change in the ratio of the semi-axes of a triple system that 
underwent adiabatic mass-loss is given by: 

Q2,f_Q2,i/' mo+mi+m 2 \ / mpj+mA 
fli.f fli.i \ OT o, f+mi +m 2 J \ mo + m\ j ' 
where a 2y \ and a 2y f are the initial and final outer semi-major 
axes, respectively, and a 2f i and a 2 j are the initial and final 
inner semi-major axes, respectively. From Equation dTTl one 
sees that the outer to inner semi-major axis ratio decreases for 
our fiducial system during adiabatic mass loss by a 2< f/a\ t f ~ 
0.8 a 2 , i/a\, i- Then with Equation (O we see that the fractional 
increase 6 in e oct due to adiabatic mass loss of the primary is 
given by: 



Eoct. f 



/ mo, 



i-mi 



\ mo — m\ 



moj+mi+m 2 \ I mo + m\ 
m.Q+mi + m 2 J \moj+mi 



(12) 



6 While e oct , ; and e ocli f will vary between systems in ^depending on ei, 
the fractional increase of e ct will, however, be the same because the mass 
loss is adiabatic so £2 drops out. 



where e octi ; and e octi f are the initial and final values of e oct , 
respectively. The absolute value is due to the ambiguity as to 
which of the binary stars will be the most massive in the final 
system. After mass loss our fiducial triple system is in the 
near-test-particle limit and e oct increases by a factor of 23.3, 
from 0.002 to 0.047. These factors enable the eccentric Kozai 
mechanism, the system then flips, and as cos i passes through 
0, (1 — «i) approaches 10 -4 . Both systems would be affected 
by tidal interactions during their first flip as r pel j < R . This 
example illustrates the essence of the MIEK mechanism. 

4. INCLINATION AND ECCENTRICITY DEPENDENCE 

The importance of the MIEK mechanism for a full popula- 
tion of triple systems will depend on the adopted mass func- 
tion of the stars, the eccentricity and semi-major axis distri- 
butions of the inner binary and outer tertiary, as well as the 
mutual inclination distribution. We save a full exploration of 
parameter space for a future paper, but we explore the effects 
of inclination and eccentricity on the dynamics of the exam- 
ple system from §f3] This system is an interesting test case be- 
cause all the masses are approximately equal, so the octupole- 
order terms mostly become important during and after mass 
loss. 

4.1. Parameters Space Exploration 

To explore the effects of inclination and eccentricity, we in- 
tegrate ~ 10 s triple systems over an equally spaced grid with 
41, 10, 20, and 4 values of cos i, e\, e 2 , and g\, respectively, 
while keeping the initial masses and initial semi-major axes 
from our example system. Because the outcomes of triple 
systems effected by MIEK may differ depending on the ini- 
tial phase angles, we integrate three triple systems at each grid 
point in cos i, e\, e 2 , and giwith randomized initial phase an- 
gles. 

In order to determine if the inner binary of each system 
becomes tidally affected or collides by our ad-hoc tidal cri- 
terion (Section I2.21 >. we record the minimum inner separa- 
tion (r m ; n ) in each of the four evolutionary phases described 
in ^2] We integrate these triple systems assuming t m \ = 10 4 
years and posit that tides will affect the inner binary if r m j n is 
smaller than some ad-hoc tidal radius (r t id e ). We define this 
ad-hoc tidal radius differently for the different evolutionary 
phases of the primary, in an attempt to capture the different 
strength of tidal forces due to the varying physical size and 
mass of the primary. Before the primary evolves off the MS 
we take r t ;d e = 3 /? ms , o, where /? mSj o in the MS radii of the 
primary. When the primary is a giant and during the mass- 
loss phase, we take r t id e = 1AU. Finally, when the primary 
becomes a WD we take r t jd e = 3 R mSy \ , where R ms \ in the 
MS radii of the secondar y in the inner bina ry. We determine 
^ mSj o and R ms ^ i from the lTorres et al.l(l2010l) sample by taking 
the median of the observed radii for stars with similar masses 
(|logM-logM'| <0.1). 

Additionally, we would like to determine when the eccen- 
tric Kozai mechanism is operating. Because orbital flips (i.e., 
a change of sign in cos z) only occur during the eccentric 
Kozai mechanism, we monitor each triple systems for flips 
during each of the aforementioned evolutionary phases. We 
then use the two above criteria to loosely define systems that 
are strongly affected by the MIEK mechanism, as systems that 
come into tidal contact and exhibit a flip after the onset of 
mass loss. 

4.2. Results 
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FIG. 2. — The "Mass-loss Induced Eccentric Kozai" (MIEK) Mechanism. Two example systems with the left and right panels having mass-loss timescales of 
'ml = 10 4 yr and f ml = 10 6 yr, respectively. Triples have m = 7.0 M Q , m x = 6.5 M Q , m 2 = 6 M , a x = 10 AU, a 2 = 250 AU, e { = 0.1, e 2 = 0.7, gi = 0°, g 2 = 180°, 
and cos i = 0.5. Top Panels: Evolution of ay (dotted black), a 2 (red dotted), e 2 (blue dotted), and mo (light blue dashed). Middle Panels: Time evolution of cos ;'. 
Bottom Panels: Time evolution of e\. Note the qualitative difference between the normal Kozai-Lidov mechanism (e\ — > 0.7) before the onset of mass loss at 3 
Myr and the eccentric Kozai mechanism (e\ 1) after mass loss when r per j < Rq (denoted by the red dashed line) and flips occur quasi-periodically. We caution 
that our integrations after the first eccentricity maximum are for purposes of illustration only since we do not include the effects of tides. 



The results of the parameter space search in inclination and 
eccentricity are presented in Table Q] Figure [3] and Figure 2] 
In Table Q] we break down the fraction of the triple systems 
that flip and become tidally affected or collide during the var- 
ious evolutionary stages presented in §|2] We highlight any 
asymmetries about cos io = by presenting the fractions for 
prograde and retrograde systems separately. We now describe 
each of the evolutionary phases of the primary in turn. 



When the primary is on the MS ("Primary MS" in Table 
[l]), 13% of triple systems become tidally affected or collide 
('"min < r tide)- We show the initial inclinations and eccentrici- 
ties of these triple systems in the left panel of Figure [3] Each 
contourplot displays the parameter space projected onto a two 
parameter plane. Each histogram marginalizes over all param- 
eters except one, and then displays the fraction of systems in 
each bin that become tidally affected or collide. From the left 
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or collided after the primary evolved off the main sequence but before the onset of significant mass-loss, r pcl j < 1AU (right). In the histograms, black shows the 
fraction of systems in each bin that are tidally affected or collide and red shows the fraction of systems that also exhibited an orbital flip. Contours present the 
fraction of systems at each grid point that would be tidally affected by our criteria. The contours are in steps of 20% starting at 10 % (black) Red dots show the 
grid points where systems where integrated. The teal line on the contours show where cos = 0. Note that there is a pronounced asymmetry between prograde 
and retrograde triples. Figure discussed in ij4] 

with e\ r max given by Equation (fl~|i and z'o the initial eccentric - 
Then, setting the periastron distance equal to r t ^ e , 



TABLE 1 
Triple Outcomes 



ity. 





Contact 


Flip 


Contact & Flip 


Unbound 


Combined 










Primary MS 


0.13 


0.064 


0.060 


0.11 


Primary Giant 


0.28 


0.000041 


0.000030 


0.00037 


During Mass-Loss 


0.00016 


0.0015 


0.000030 


0.0027 


After Mass-Loss 


0.022 


0.021 


0.019 


0.011 


Prograde 










Primary MS 


0.039 


0.0076 


0.0059 


0.15 


Primary Giant 


0.30 


0.000063 


0.000042 


0.00042 


During Mass-Loss 


0.00013 


0.0020 


0.000021 


0.0037 


After Mass-Loss 


0.030 


0.029 


0.028 


0.0083 


Retrograde 










Primary MS 


0.22 


0.12 


0.11 


0.082 


Primary Giant 


0.26 


0.000020 


0.000020 


0.00032 


During Mass-Loss 


0.00020 


0.00097 


0.000040 


0.0017 


After Mass-Loss 


0.015 


0.013 


0.010 


0.014 



NOTE. — Fractions of triple systems that become tidally affected, exhibit a flip, or 
are unbound at various evolutionary stages without being tidally affected or unbound in 
a previous evolutionary stage. Prograde and retrograde divide the sample in half where 
combined are the fractions of the whole sample. The rows Primary MS, Primary Giant, 
During Mass-Loss and After Mass-Loss present the fractions of triple systems with the 
primary star is on the main sequence, a giant that has evolved off the main sequence, 
an evolved star undergoing an 10 4 year phase of significant mass-loss phase and a WD, 
respectively. The columns are as follows. Contact — Systems that satisfy our ad-hoc 
tidal criteria. Flip — Systems whose cos i changes sign from its initial value. Contact 
& Flip — Systems common to the previous two columns. Unbound — Systems that 
become unbound without being tidally affected prior. Table discussed in ^4] 



panel of Figure [3] we see that when the primary is on the MS 
the majority of systems which become tidally affected or col- 
lide are initially highly inclined. These systems undergo nor- 
mal Kozai-Lidov oscillations, which bring the inner binary to 
tidal contact with a maximum pericenter distance given by: 



' pen 



: <3l(l-<?l,max) = ai 



1- 



1 J 2 ■ 
1 COS In 

3 ° 



1/2' 



(13) 



' pen ' 



■ Aide = 3fl ms , ~ 12R - 0.05AU - 0.005 a b (14) 

we see that the required critical initial inclination to interact 
tidally when the primary is on the MS is 



cos/o,crit ~ 0.08, 



(15) 



which is in good agreement with the prograde orbits in the 
left panel of Figure [3] Note, however, that the distribution is 
asymmetric about cos io = 0, with 22% of the retrograde sys- 
tems becoming tidally affected or colliding while only 3.9% 
of the prograde systems do. This asymmetry is not captured 
by Equation (fTJ and Equation ( fT3l ), because of the approxi- 
mations under which Equation (d} is derived. A more com- 
plete treatment of the quadrupole-order Hamiltonian shows 
that ei.max does not occur at cos i n = 0, but instead at retrogr ade 
mutual inclination (see Fig. 1 of Miller & Hamilton 2002). 

Importantly, we find that 6.0% of the triple systems, with 
cos /o ~ 0, become tidally affected and exhibit an orbital 
flip, which is characteristic of the eccentric Kozai mecha- 
nism. These systems comprise almost half of the total num- 
ber of systems that become tidally affected or collide when 
the primary is on the MS, showing that the eccentric Kozai 
mechanism is also important. The implications that the ec- 
centric Kozai has for previous studies involving the Kozai- 
Lidov mechanism and tidal friction is discussed further in §|5l 
but here we note that the systems that tidally interact or col- 
lide on the MS are an a logou s to the systems investigated by 
Fabry cky & Tremaine (2007), except that these authors only 
included quadrupole-order terms in the secular Hamiltonian, 
and thus may have missed the systems that undergo the ec- 
centric Kozai mec hanism and execute a "flip" with e\ — >■ 1 
(iNaozet alJl2^Tilbh . 

Additionally, 11% of triple systems become unbound 
when the primary is on the MS. The systems that are un- 
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FIG. 4. — Initial parameters for triple systems that were tidally affected or collide after the end of mass loss, r pcr j < 3R ms [ . In the histograms, black shows the 
fraction of systems in each bin that are tidally affected or collide and red shows the fraction of systems that also exhibited an orbital flip. Contours present the 
fraction of systems at each grid point that would be tidally affected by our criteria. The contours are in steps of 20% starting at 10% (black). Red dots show the 
grid points where systems where integrated. The teal line on the contours show where cos = 0. Note there is a pronounced asymmetry between prograde and 
retrograde triples. Figure discussed in 



bound have large ei and do not satisfy the standard, empir- 
ically derived, stability cri terion for triple systems given by 
( Mardling ^&Aarsethll^OOll) : 



a 2 



>Cf 



1 + - 



/)?2 



l+e 2 



OTo + OTi I (1 — ei) 



0.4 



(16) 



where C = 2.8 and f = 1 - —i is an ad-hoc inclination term. 
These systems are present in our study because we uniformly 
sample the eccentricities without regard for the initial dynam- 
ical stability of the system. We see from Equation dT6b and 
our initial assumed semi-major axes, that if cos i = then 
<?2 < 0.83 must hold for the system to satisfy the stability 
criterion. It is interesting to note that in the cos ;'o vs. <?2, o 
contour plot in the left panel of Figure [3] many triple systems 
which are unstable given the stability criteria become tidally 
affected or collide before they are unbound. It is possible that 
tidal friction will quickly shrink a\ of these systems, causing 
them to become stable as the ratio a^jax increases. 
After the primary has evolved off the MS but before signif- 



icant mass-loss has begun ("Primary Giant" in Table [TJ, 28% 
of triple systems become tidally affected or collide (r m ; n < 
1AU) for the first time. We show the initial inclinations and 
eccentricities of these triple systems in the right panel of Fig- 
ure [3] Such a large fraction of the triple systems become 
tidally affected or collide during this stage of evolution be- 
cause the physical size of the primary dramatically increases, 
so we impose an ad-hoc increase to r t jd e = 1AU. Then from 
Equation ( fT3l ) we find that 



cos /o, crit ~ 0.34, 



(17) 



encompassing about a third of our parameter space. Addition- 
ally, we see that only a small fraction of the systems exhibit 
the orbital flip characteristic of the eccentric Kozai mecha- 
nism when the primary is a giant, showing that the eccen- 
tric Kozai mechanism is not important during this phase of 
the primary's evolution. Lastly, only 0.037% of the systems 
are unbound during this evolutionary phase, emphasizing that 
initially dynamically unstable systems are either quickly un- 
bound, when the primary is still on the MS, collide, or become 
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tidally affected. 

During the f m i = 10 4 year mass-loss phase ("During Mass- 
Loss" in TableQ} only a small fraction, 0.016%, of the triple 
systems become tidally affected or collide for the first time. 
This is because previous stages of evolution have removed the 
majority of triple systems that would tidally interact or col- 
lide with r t ide = 1 AU, and because t m \ < t&, so the systems do 
not have time to complete a Kozai cycle during the primary's 
mass loss episode. We find that the fraction of systems which 
are unbound increases, to 0.27%, from the previous phase of 
evolution. More systems become unstable because, as shown 
by Equation ( fTTT ). the outer to inner semi-major axis ratio de- 
creases during mass loss. Then, as a result of Equation ([TBI , 
this decrease in the semi-major axis ratio causes more triple 
systems to be unstable. 

Finally, after the end of mass-loss ("After Mass-Loss" in 
TableUJ we see that 2.2% of triple systems spanning a broad 
range in orbital parameters become tidally affected or collide 
for the first time. We note that the lion's share, 1 .9%, of these 
systems exhibit an orbital flip characteristic of the eccentric 
Kozai mechanism and thus the MIEK mechanism because it 
is induced by mass loss of the primary. We now discuss the 
distribution of tidally affected or colliding systems after mass 
loss for the inclination and eccentricities individually. 

The distribution of these tidally affected or colliding sys- 
tems is asymmetric about cos ;'o = 0, with 3.0% of pro- 
grade systems and only 1.5% of retrograde systems becom- 
ing tidally affect of colliding. At least some of this asym- 
metry arises because, like previous stages of evolution of the 
primary, systems with cos io near are strongly preferred. 
However, earlier stages of evolution have removed a dispro- 
portionate fraction of the prograde systems, 33%, compared 
to the fraction of retrograde systems, 49%. This causes the re- 
verse asymmetry, more prograde than retrograde systems, in 
the distribution of systems which become tidally affected or 
collide after the end of mass loss. 

The distribution in <?2, o for the systems that become tidally 
affected or collide after the end of mass loss shows a strong 
preference for larger values, peaking at ei. o = 0.75. This pref- 
erence is easily understood when considering Equation (f2), 
which also shows a strong dependence on <?2, o- Thus, the 
octupole order terms and consequently the eccentric Kozai 
mechanism are expected to become more important with 
larger <?2, o- Even though there is a strong preference for large 
e2, o, the MIEK mechanism continues to operate for systems 
with modest <?2, o ~ 0.4. We also note that the <J2, o distribution 
is truncated at very large values, which arises because most 
of these systems became tidally affected, collided, or were 
unbound when the primary was on the MS. 

Lastly, the distribution of e\ t o for the systems that become 
tidally affected or collide shows a preference for both large 
and small values. The minimum in the distribution occurs at 
ei,o ~ 0.5, with ~ 3 times fewer systems becoming tidally 
affected or colliding then at e\, o = 0. The preference for large 
values of ei,o is easily understood because these systems do 
not need to increase their inner eccentricity as much to collide 
or come into tidal contact. The preference for smaller values 
of ei.o, however, is not easily understood. We also note that 
the majority of systems with o = 0.9 become tidally affected 
or collide in earlier stages of the primary's evolution, so these 
systems are absent in the current ej.o distribution. 



We now briefly evaluate the sensitivity of the results pre- 
sented in 34.21 on the ad-hoc tidal criterion chosen (Section 

First, we evaluate the dependence of the fraction of systems 
that become tidally affected or collide on our choice of tidal 
criterion. We both decreased and increased the tidal criterion 
from r tide = 3 R* to r tide = 2.5 /?», r tide = 3.5 /?», r tide = 5 R*, 
and r t ide = 10.0 R* = R ms . o when the primary is on the MS 
and /?„ = ^ ms> i when the primary is a WD. The fraction of 
systems that become tidally affected or collide after the end 
of mass-loss with these new tidal criteria are 0.022, 0.023, 
0.023, and 0.025, respectively. When these fractions are com- 
pared to that presented in Table Q] we see that the fraction of 
systems affected by the MIEK mechanism is relatively insen- 
sitive to the exact choice of tidal criterion. This insensitivity 
results because an increase in the tidal radius on the MS af- 
fects systems that will already be affected by the large tidal 
radius when the primary becomes a Giant, and does not af- 
fect regions of parameter space shown in Figure |4] where the 
MIEK mechanism operates. In fact, increasing the tidal ra- 
dius marginally increases the fraction of systems affected by 
MIEK. 

Additionally, we emphasize that by eliminating all systems 
that collide or were tidally affected at a prior stage of evolu- 
tion, 1 .9% is a conservative estimate of the fraction of triple 
systems that will undergo the MIEK mechanism. This is be- 
cause many of the systems we remove from the sample as the 
evolutionary stages progress may still be viable candidates for 
the MIEK mechanism. One possibility is that tides on the gi- 
ant primary only detune the normal Kozai-Lidov mechanism, 
and do not allow e\ to increase so dramatically. Then, our 
removal of these systems from the sample would unrealisti- 
cally decrease the number of systems that experience MIEK. 
These systems would survive in our sample to the stage after 
mass loss is complete. Another possibility is that the com- 
bination of the Kozai-Lidov mechanism and tides on a giant 
will work to shrink a\ to a few AU. These systems will then 
have 02/01 ~ 100, but they can still be affected by the eccen- 
tric Kozai mechanism after mass loss, which will shrink this 
semi-major axis ratio according to Equation (fTTt . 

Evaluating these effects would involve a realistic calcula- 
tion of the 3-body dynamics with tides since cos i will evolve 
secula rly during the tidal interaction (Fabr vckv & Tremaind 
12007b . We save this for a future work. Here, to roughly quan- 
tify the number of systems that could potentially be affected 
by the MIEK mechanism, we changed the tidal criterion when 
the primary star is a giant and during mass loss to be the same 
as it was when the primary star was on the MS, r t j de = 3/? mSi o. 
Under this new tidal criterion, the analogue of Figure |4]is pre- 
sented in Figure [5] As expected, we see a sharp decrease, 
from 28% to 0.01 1%, in the fraction of triple systems that are 
tidally affected or collide when the primary is a giant. There 
are then many additional systems available after mass loss to 
increase the fraction of systems that are tidally affected or col- 
lide after the primary becomes a WD. The total change in the 
fraction is from 2.2% to 8.7%. Similarly, there is a significant 
increase in the fraction of MIEK systems, systems that col- 
lide or become tidally affected and flip the sign of cos i, from 
1.9% to 6.7%. For the example system considered, and for a 
uniform distribution in eccentricities and mutual inclination, 
we consider 1.9% and 6.7% to be lower and upper limits, re- 
spectively, to the fraction of all systems that undergo MIEK. 



4.3. Sensitivity to r t id e 



5. DISCUSSION 
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Fig. 5. — Same as Figure [4] except assuming the primary's tidal radius 
remains unchanged as it evolves until it becomes a WD. This is taken to be 
an upper bound on the fraction of systems that become tidally affected or 
collide after the end of mass loss. Figure discussed in !j4] 

It is interesting to speculate on the outcome of triple sys- 
tems that become tidally affected at each of the various points 
in their evolution. In what follows, we discuss the possible 
outcomes for tidally interacting triple systems during the MS, 
giant, and WD evolutionary phases of the primary in §15.11 
jj5.2| and 35.31 respectively. We then discuss the MIEK mech- 
anism for massive star triples during NS formation of the pri- 
mary in {s!5.4l Lastly, we discuss extensions of the current 
study in ^531 

5.1. Main Sequence Phase 

We have shown in our fiducial parameter search, that sys- 
tems that interact tidally while the primary is on the MS un- 
dergo both normal Kozai-Lidov oscillations and the eccentric 
Kozai mechanism in approximately equal proportion. Those 
systems that undergo normal Kozai-L idov oscillations are 
analog ous to the systems considered by Fabr ycky & Tremainel 
(2007), who calculated the long-term evolution of a dis- 
tribution of triple systems using the quadrupo l e-orde r 
secular equations of lEggleton & Kiseleva-Eggletonl (1200 ll) . 
iFabrvckv & Tremainel (12007b found that the combination of 
normal Kozai-Lidov cycles and tidal friction produces a sub- 
stantial population of close binarie s, and that their calcu - 
lations matched the observations of iTokovinin et al.l d2006l) . 
which show that essentially all close spectroscopic binaries 
are actually triple systems. Given these results, we anticipate 
that our systems with cos z'o ~ will also tidally circularize 
with small a\ and P\, with the caveat that the MS lifetime of 
the primary we cons ider here is much shorter than for the sys- 
tems considered by Fabrycky & Tremaine (2007), and thus a 
fraction of our systems may not circularize before WD forma- 
tion. 

Importantly, though, we find that half of the systems that 
tidally interact or collide on the MS do so as a result of 
the ec ce ntric Kozai mecha nism, which, as iLithwick & Naoj 
d20ll . iKatz et alJ ( I20T1 . and iNaozetalJ ( 12011 alibi) have 
shown c annot be captured by the qua drupole-order calcula- 
tions of IFabrvckv & Tremaine] (120071) . What is the fate of 
these systems? Unfortunately, no study has yet been per- 



formed to investigate the eccentric Kozai mechanism with the 
inclus ion of tida l diss ipation for triple stellar systems. How- 
ever, iNaoz et all (1201 lal) demonstrate with an octupole-order 
calculation that the eccentric Kozai mechanism and tidal fric- 
tion can possibly explain the occurrence of retrograde hot 
Jupiters. They show that the extremely high eccentricities ob- 
tained by the inner binary when cos i flips signs can lead to 
a rapid capture of a planet into a short period retrograde orbit 
("Kozai capture."). We suppose a similar mechanism oper- 
ates for stellar triple systems, and that those that are tidally 
affected on the MS due to the eccentric Kozai mechanism 
rapidly dissipate energy from the inner binary and circularize 
with a semi-major axi s of a\ « 2r per j ~ 2r t {& e . If so, then the 



IFabrvckv & Tremainel (120071) study may be missing of order 
half of the systems that become close binaries. 

The action of the eccentric Kozai mechanism for MS bi- 
naries may also affect the distribution of binary eccentricities 
as a function of inner period as th e population is circularized 
(ISocrates et al.ll201 ll Dong et al.ll2012l) . 

An example observed (probable) intermediate mass triple 
star system that has likely already tidally circularized is 
MT429 in the Cygnus OB2 association, where the inner Pi ~ 
3 d B3V+B6V binary is orbited by a B0V tertiary at 138 AU 
dKiminki et alJ [2012). An example of an even more mas- 
sive system that might give rise to a NS-MS triple system as 
discussion in Section [5~4l is HD 150136. It consists of 03- 
3.5V+05.5-6V (- 64 + 40 M Q ) inner binary with Pi ~ 2.7 d 
and an Q6.5-7V ( ~ 35 M ) tertiary with P 2 ~ 3000-5000d 
dMahv et al.ll2012h . 

Finally, after the inner binary orbit has shrunk due to either 
normal Kozai-Lidov oscillations or the eccentric Kozai mech- 
anism, the binary evolves like other close binaries, potentially 
leadin g to the progenitors o f single degenerate (SD) super- 
novae (Whelan & Iben ll 19731). double degener ate (DD) super- 
novae J^bbinjiniM^ barium stars 
(BA; iMcClure et alTT l980). and asymmetric planetary nebu- 
lae (PNT^Monlliasl - 

5.2. Giant Phase 

We find that ~ 30 % of the triple systems we consider in- 
teract tidally or collide as the primary evolves off the MS, but 
before mass loss, mostly as a result of normal Kozai-Lidov 
oscillations, which bring the pericenter of the secondary or- 
bit inside ~ AU (see right panel, Fig. [3] and Equation [TTb . 
Since the post-MS phase is short (particularly for the inter- 
mediate and high mass stars we consider here), and since the 
efficiency of tidal friction is uncertain, it is unclear if the sys- 
tem will circularize before mass loss, or whether the Kozai- 
Lidov oscillation will simply be de-tuned such that pericen- 
ter ~ AU for the duration of the giant phase. The interaction 
of Kozai-Lidov oscillations as the primary evolves to the gi- 
ant phase may also initiate mass transfer and complex sec- 
ula r dynamics, perhaps ana logous to the system considered 
bv lProdan & Murravl (120121) . Another possibility is that this 
phase initiates common envelope evolution, perhaps circular- 
izing the inner orbit at smaller semi-major axis. If the inner bi- 
nary maintains significant eccentricity, or if ei is pumped after 
the giant phase and during mass loss via normal Kozai-Lidov 
oscillations or MIEK, such systems ma y form off-center 



or bi-polar planetary nebulae (PNe) (e.g., Soker et ail 119981 : 
ISoker & Rappaportl l2000t IWitt et all |2009|). On the other 
hand, if the inner binary is circula rized, it might also par- 
ticipate in t he shaping of PNe (e.g., [Mastrodemos & Morris 
I1998LI1999I) . An example of some of these mechanisms in ac- 
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FIG. 6. — MIEK for massive stars, similar to Figure|2] but here mo = 11 Mq , 
mi = 10M Q , m 2 = 1OM , e\ =0.7, e 2 =0, cosi = 0.25, gi = 0°,g 2 = 270°, 
a\ = 15 AU, and a 2 = 300 AU. At ~ 7 Myr, mo decreases to 1.4 Mq, cos/ 
flips, and the system attains 1— e\ ~ 10~ 3 immediately. 

tion may be ring planetary nebula SuWt 2 (Exter et al. 2010). 
Since a very large fraction of triple stars will be first affected 
by Kozai-Lidov cycles when the MS star evolves, this phase 
should be investigated in detail. 

5.3. After Mass Loss & MIEK 

We find that ~ 2 % of our triple systems undergo the MIEK 
mechanism after WD formation (see Fig. |2). We again em- 



phasize that by removing all systems from our sample that 
were tidally affected during the giant phase, the fraction of 
triples that go through MIEK presented in Table Q] ("After 
Mass-Loss") is conservative since many of these systems may 
be available to go through MIEK. The reason for this is that 
as the primary enters its post-MS evolution and begins to in- 
teract tidally with the secondary, the orbits will be affected 
by pr ocesses that we do not model (e.g., iProdan & Murravl 
120121) . We attempt to provide an upper limit on the fraction 
of triple systems that could be tidally affected or c ollid e after 
mass-loss due to the MIEK mechanism in Section l4~3l bv as- 
suming a small value for r n d e during the giant phase (compare 
Figs. 2] and El. This increases the fraction of triple systems 
that undergo MIEK from 1.9% to 6.7%. A calculation of the 
3-body dynamics with tides appropriate to a MS star interact- 
ing with a red giant is clearly needed to address this issue in 
more detail and to address the issues raised in Section 15721 

Systems that undergo the MIEK mechanism also have many 
possible outcomes. One possibility is that the extreme eccen- 
tricities obtained during a "flip" might lead to physical colli- 
sions (Baoz et al. in prep.). This may occur if, in a single orbit 
of the inner binary, r per i changes from r per j > r t id e to r per i < 
This is possible because the angular momentum of the inner 
binary at such large e\ is small. Thus, the change in angu- 
lar momentum required for r per i ~ is also small. Another 
possible out come for the system is a tidal Kozai capture as 
presented in iNaoz et al.l d201 lal) . The inner semi-major axis 
«i will strongly decrease because of tides, leaving a close MS- 
WD binary. When the secondary subsequently evolves off the 
MS, the WD will either accrete or go through common enve- 
lope evolution. Both of these scenarios could lead to single- 
or double-degenerate (SD or DD) SNe la, cataclysmic vari- 
ables, or other types of WD accretors such as AM CVn stars 
(Warner] 1 19951). whi ch may produce faint ".la" supernovae 
Bildsten et afl|2007l) . 

One interesting consequence of this evolutionary scheme 
for producing close WD-MS binaries is that it skips the initial 
common envelope evolutionary phase usually required in bi- 
nary evolution to produce a close WD-MS pair. In the MIEK 
mechanism the close binary is produced without common en- 
velope, and thus the WD produced should have a mass ap- 
propriate to a sing le star with the mass of the primary (e.g., 
iKalirai et al.l 120081) . This stands in contrast to the conven- 
tional picture for the production of WD-MS and WD-WD bi- 
naries relevant for SD and DD la progenitors, where the WD 
growth is truncated by common envelope interaction with the 
close secondary and the mass of the resulting WD is smaller 
than one would naiv ely estimate from the single-star initial- 
final mass relation dlben & Tutukovlll984l) . For SD SN la 
models, this means that less mass needs to be accreted onto 
the WD from the secondary, leading to a shorter delay time 
between production of the system and explosion. For DD 
SN models, if the model presented in Pakmor e t al.l (1201 2b 
is accurate, then the luminosity of the SN resulting from the 
merger of a WD-WD binary is primarily determined by the 
mass of the primary. Thus, all else being equal, if a WD-WD 
binary results from a MIEK triple system then its primary WD 
mass should be larger and thus result in a more luminous DD 
SN than the analogous close binary system that is produced 
through normal common envelope evolution. 

5.4. Massive Star Triples & Neutron Star Formation 

The MIEK mechanism may also play an important role in 
triple systems after NS formation, particularly since the NS 
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mass, ~ 1.4 M©, is so much less than its progenitor massive 
star, leading to a large increase in e oc t at the time of the associ- 
ated supernova (Equation |2). However, several effects differ 
in the NS case versus the case of WD formation considered 
throughout this paper. First, NSs receive a "kick" at birth that 
may unbind the inner binary or the outer tertiary, depending 
on the masses of t he constituents, their eccentr icities, and their 
semi-major axes (lHillslll983t lKalogeral ll996). Second, mass 
loss on the MS can affect the mass ratios as a function of time 
(e.g jEkstrom et ail201 ll) . Third, even in the absence of a NS 
kick, the large and instantaneous mass loss in the supernova 
explosion from the primary provides a kick to the center of 
mass of the inner binary which may unbind the tertiary, and 
shift the mutual inclination. Even so, we expect a popula- 
tion of NS-MS binaries to survive with massive tertiaries and 
large e 2 due to the aforementioned kicks. Thus, e oct will be 
large for the surviving triple systems which will cause many 
to exhibit the MIEK mechanism. This is a qualitatively new 
way to produce close NS-MS binaries. 

As a test, we integrated the orbit of a system with mo = 
1 1 M Q , mi = 1OM , m 2 = 10M Q , e x = 0.7, e 2 = 0, cos/ = 0.25, 
gl = 0°, g 2 = 270°, fli = 15 AU, and a 2 = 300 AU. At a time 
~ 7 Myr after the start of the simulation we instantaneously 
decrease the mass of the primary to 1 .4M Q and apply no kick 
to the NS. Figure|6]shows the resulting evolution. The system 
flips immediately after the supernova and evolves to an inner 
eccentricity of 1 —e\ > 10~ 3 , bringing the NS to tidal contact 
with the secondary. Note the large increase in e 2 from ~ 
to ~ 0.3 immediately after the SN. This calculation is meant 
only to be illustrative since it does not include a NS kick. A 
full population study, with kicks and tides is clearly needed to 
assess the viability of this mechanism for forming close NS- 
MS binaries. 

A further aspect of this type of evolution is that the ratio of 
a 2 ja i can decrease during the strong MS mass loss and the 
supernova explosion, which can cause the triple to become 
dynamically unstable and eventually disru pt, or potentially 
initiat e collisions (the TEDI mechanism of Perets & Kratter 
l2012h . 

5.5. Further Studies 

In this study we have explored a single set of mo and m\, 
however, different combinations of masses will affect the frac- 
tion of systems which become tidally affected or collide dur- 
ing each of the evolutionary phases. For example, if mo ~ m\ 
the system will have little or no time after the primary be- 
comes a compact object before the secondary enters its giant 
phase. This will reduce the fraction of systems that will flip 
and become tidally affected or collide due to the MIEK mech- 
anism. However, if the difference between the masses of the 
inner binary members is significantly greater, then a larger 
fraction of systems will flip before the primary goes through 
mass-loss. Additionally, the difference in the masses after the 
primary has become a compact object will be smaller, further 
decreasing the fraction of systems that will be affected by the 
MIEK mechanism. Determining the optimal balance between 
these two effects is beyond the scope of the current study. 

In a future study we will present a more general study of 
mass loss in a full population of triples, quantifying the num- 
ber of systems affected by Kozai-Lidov oscillations and the 
eccentric Kozai mechanism, an accounting of those that are 
unbound and those that are tidally affected or collide at each 
evolutionary stage, and those that undergo the MIEK mecha- 
nism. 



Any study which evolves a full population of triple sys- 
tems with distributions in a\ and a 2 will have a signifi- 
cant fraction of s ystem s that will be affected by the TEDI. 
iPerets & Kratteri (1201 2|) demonstrate that systems that start 
with a small 02/01 ratio become unstable during mass-loss 
and evolve chaotically, leading to close encounters, collisions, 
and exchanges b e tween the stellar components. Additionally, 
IPerets & Kratteri d2012l) also raise the issue of Kozai cycles 
and stellar evolution as important. Future studies will be 
needed to quantitatively evaluate the importance of both the 
TEDI and the MIEK mechanism. 



6. CONCLUSION 

We have presented a preliminary investigation of stellar 
evolution and mass loss in triple systems and demonstrated 
the MIEK mechanism, a novel channel through which triple 
systems produce compact object - MS binaries that interact 
tidally or collide only after the primary has evolved off the 
MS and become a compact object. 

For a broad range of parameters, hierarchical triple star sy s- 
tems are unaffected by Kozai-Lidov oscillations until the pri- 
mary in the central binary evolves off the MS and begins mass 
loss. Subsequently, the primary becomes a WD or a NS, and 
may then be much less massive than the other components in 
the ternary, enabling the "eccentric Kozai mechanism". In this 
near-test-particle limit, the mutual inclination between the in- 
ner and outer binary can flip signs, driving the inner binary 
to very high eccentricity and tidal contact or collision. In this 
study, we define tidal contact by an ad-hoc minimal separation 
between the member of the inner binary, below which tidal af- 
fects are deemed important (Section [2.2b . Even distant bina- 
ries with initial semi-major axes larger then tens of AU can be 
strongly affected. We consider an example triple system with 
masses 7, 6.5 and 6 M Q as proof -of-principle (see Fig. |2}, 
and explore the MIEK mechanism's dependence on the initial 
eccentricities and inclination for this system. For a flat distri- 
bution of eccentricities and cos/, we find that 13%, 28%, and 
2.2% of systems are tidally affected or collide before the pri- 
mary evolves off the MS, after the primary becomes a giant, 
and after the primary becomes a WD, respectively. 

On the MS, we find that roughly half of our systems inter- 
act tidally or collide as a result of the eccentric Kozai mecha- 
nism, and hence their dynamics would not have been captured 
by quadrupole-order secular calculations. In the giant phase, 
most of the systems are brought to tidal contact by normal 
Kozai-Lidov oscillations and the fact that pericenter passage 
occurs at semi-major axes less than ~ AU even for very mod- 
est inclinations. The large fraction of such systems motivates 
a detailed study of the dynamics and tidal interaction of triple 
systems with a giant primary. Finally, the MIEK mechanism 
dominates the last stage, causing the WD produced by the 
initial primary to come to tidal contact or collide with the ini- 
tial secondary which reduces a\ without a common envelope 
phase. As a last application, we showed that some massive 
star triple systems will also undergo the MIEK mechanism af- 
ter the primary's supernova explosion and the associated rapid 
mass loss (Figure |6}. 

We save a detailed study of the parameter space of triples 
affected by MIEK for a future paper, but here note that for a 
thermal distribution of eccentricities, as is commonly adopted 
for studies of binary stars, we expect more systems to be af- 
fected by MIEK since, in general, e oct increases as e 2 increases 
(Equation|2l). 
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